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We calculate the dominant Standard Model contributions to the time-dependent CP asym- 
metry in B° — > i^*°7, which is 0(l/m&) in QCD factorisation. We find that, including all 
^ ■ relevant hadronic effects, in particular from soft gluons, the asymmetry S is very small, 
c3 ' S 1 = —0.022 ± 0. 015+q 01 , and smaller than suggested recently from dimensional arguments 
in a 1 /nib expansion. Our result implies that any significant deviation of the asymme- 
try from zero, and in particular a confirmation of the current experimental central value, 
5hfag = —0.28 ± 0.26, would constitute a clean signal for new physics. 
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1 Introduction 



The radiative decay b — > 57 has been extensively studied as a probe of both the flavour 
structure of the Standard Model (SM) and new physics beyond the SM (see Ref. [1] for 
a review). While the vast majority of studies has focused on the prediction of the decay 
rate for exclusive and both spectra and decay rate for inclusive b — > 57 decays, there is one 
rather peculiar feature of this process which has attracted far less attention, namely that, 
in the SM, the emitted photon is predominantly left-handed in b, and right-handed in b 
decays. This is due to the fact that, in the language of effective field theories, the dominant 
contribution is from the chiral-odd dipole operator SL^a^bn^y As only left-handed 
quarks participate in the weak interaction, an effective operator of this type necessitates, 
in the SM, a helicity flip on one of the external quark lines, which results in a factor rrib 
(and a left-handed photon) for b^ — > sl7l and a factor m s (and a right-handed photon) for 
fex, — ► s_r7_r. Hence, the emission of right-handed photons is suppressed by roughly a factor 
m s /m} } . This suppression can easily be alleviated in a large number of new physics scenarios 
where the helicity flip occurs on an internal line, resulting in a factor mi/m b instead 
of m s /m h . A prime example are left-right symmetric models [2], whose impact on the 
photon polarisation was discussed in Ref. [3]. These models also come in a supersymmetric 
version whose effect on b — > S7 was investigated in Ref. [4]. Supersymmetry with no left- 
right symmetry can also provide large contributions to b — > S7^j, see Ref. [5] for recent 
studies. Other potential sources of large effects which have been studied are warped extra 
dimensions [6] or anomalous right-handed top couplings [7]. Unless the amplitude for 
b — > S7# is of the same order as the SM prediction for b — > sj L , or the enhancement of 
b — > S7r goes along with a suppression of b — > s^l, the impact on the branching ratio 
is small, as the two helicity amplitudes add incoherently. This implies there can be a 
substantial contribution of new physics to b — > S7 escaping detection when only branching 
ratios are measured. 

Although the photon helicity is, in principle, an observable, it is very difficult to measure 
directly. It can, however, be accessed indirectly, for instance in the time-dependent CP 
asymmetry in B° — ^ K* 0/ ~f, which relies on the interference of both left and right helicity 
amplitudes and vanishes if one of them is absent. This method was first suggested in 
Ref. [3] and later discussed in more detail in Refs. [8, 9]. It is rather special in the sense 
that usually new physics modifies the SM predictions for time-dependent CP asymmetries 
by affecting the mixing phase (as in B s — > J/ip(f), see for instance Ref. [10]), introducing 
new weak phases or moderately changing the size of the decay amplitudes which, in the 
absence of precise calculational tools, makes it very hard to trace its impact. In contrast, 
the time-dependent CP asymmetry in B° — > K* 0/ ~f is very small in the SM, irrespectively 
of hadronic uncertainties, by virtue of the helicity suppression of one decay amplitude, 
and new physics enters by relieving that suppression. The smallness of the asymmetry in 
the SM, and the possibility of large effects from new physics, makes it one of the prime 
candidates for a so-called "null test" of the SM, as recently advertised in Ref. [11]. 

Other channels and methods to probe the photon helicity have been discussed in 
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Figure 1: Dominant contribution to b — ► s^/g. A second diagram with photon and gluon vertices 
exchanged is implied. 

Refs. [12, 13]. In this letter, however, we focus on the time-dependent CP asymmetry 
in B° — > K* ^. It is given by 

where and -ft^* are observed via their decay into the CP eigenstate Ks7r°. The term 
involving an interference of photons with different polarisation is S, for which the following 
experimental results are available from the B factories: 

^BaBar = -0.21 ± 0.40 (stat) ± 0.05 (syst) BaBar [14] (232 • 10 6 BB pairs), 
^Bciie = — 0.32lg;|| (stat) ± 0.05 (syst) Belle [15] (535 • 10 6 BB pairs), (2) 

with the HFAG average Shfag — —0.28 ± 0.26 [16]. While these results are compatible 
with zero at the la level, the central values of both BaBar and Belle are in agreement 
and interestingly large. A drastic reduction of the experimental uncertainty will probably 
be difficult at the LHC, but can be achieved at a Super B factory, with an anticipated 
statistical uncertainty of S of 0.07 with 10 ab" 1 of data [17] and 0.04 with 50 ab" 1 [18]. 

In order to clearly distinguish any new physics signal from the SM background, one 
needs to know the latter as precisely as possible. As discussed above, one contribution 
comes from b^ — > s_r7r, with a helicity flip on the s quark line; it generates the contribution 

S SM ' S " = - sin(2/3) ^ (2 + 0(a s )) (3) 
m b 

to the CP asymmetry, with (3 being one of the angles of the CKM unitarity triangle. At 
leading order in a s , S SM ' SR is free from hadronic uncertainties. As pointed out in Ref. [8], 
another mechanism to remove the helicity suppression of b — > is to emit an additional 
gluon. The dominant contribution to this mechanism is via a c-quark loop and is shown 
in Fig. 1. In inclusive decays this is a bremsstrahlung correction and can be calculated in 
perturbation theory [8]. In exclusive decays, on the other hand, the gluon can be either hard 
or soft. If it is hard, it attaches to the spectator quark, which induces 0(a s ) corrections to 
(3). If it is soft, it has to be interpreted as a parton in one of the external hadrons. Stated 
differently, if the gluon is soft, the amplitude involves higher Fock states of the B and K*. 
A data-driven method to distinguish this contribution from that of the dipole operator 
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SL(R)0'^bR(L) was discussed in Ref. [13] and relies on the Dalitz-plot analysis of decays of 
type B° — > ^K s + neutrals, where neutrals stands for ir°, rf\ K s , light vector mesons or 
any combination of these particles. In Ref. [8] it was shown, in the framework of soft- 
collinear effective theory (SCET), that contributions from b — > s^yg are suppressed by one 
power of rrib with respect to the left-handed photon emission, which confirms the results 
obtained in QCD factorisation [19, 20], where an explicit 0{a s ) calculation demonstrated, 
to leading order in 1/rrib, the absence of right-handed photons in B° — > K* 0/ ~f. A SCET 
analysis of the CP asymmetry in Ref. [9] estimated the size of the l/m b corrections to S SM 
induced by b — > sjg as ~ ±0.1, but is based on dimensional counting of the operators 
involved rather than a calculation of the relevant matrix elements. Another calculation, in 
perturbative QCD, gives S^q CD = —(3.5 ± 1.7)% [21], including effects mainly from hard 
gluons; the contribution of soft gluons is treated in a model-dependent way. 

The purpose of this letter is to provide a calculation of the soft-gluon contribution 
to the time- dependent CP asymmetry in B — > K*^y induced by the c-quark loop shown 
in Fig. 1. The method we use are QCD sum rules on the light-cone. It turns out that 
the relevant hadronic parameters were calculated already previously, in 1997, in Ref. [22], 
using the method of local QCD sum rules. The motivation at the time was to estimate 
long-distance corrections to the branching ratio of B — > K*^. In fact those corrections 
were first discovered for the inclusive process [23]. In this letter, we show that the same 
parameters also enter the time-dependent CP asymmetry in B — > K*^ and present a new 
calculation of their values. 



2 The CP Asymmetry 

Let us define the amplitudes of the decay of B mesons into K* and left- or right-handed 
photons in the following way: 

Al { r) = A(B° - K*° lL{R) ) , Al { r) = A(B° - K*° lL(R) ) . (4) 

Neglecting, as usual, the small width difference between B° and B°, the time-dependent 
CP asymmetry is then given by (1) with 
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Here q/p is given in terms of the B°-B° mixing matrix M 12 , in the standard convention 
for the parametrisation of the CKM matrix, by 



p V M 12 

Extending in an obvious way the notations introduced in Ref. [20] in the context of 
QCD factorisation, the decay amplitudes can be written as 

Al(r) = ^ [Ka u 7 {K* lL{R) ) + \ c a c 7 (K* lL{R) )) {K* lL(R) \Q L 7 {R) \B) 
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= ^ {\ u a u 7L(R) + X c a c 7L{R) ) (K* lL(R) \Q L{R) \B) . (6) 

In QCD factorisation, a 7 2 are of order 1 in a l/m b expansion [20], 1 

c$2 = C 7 + 0(a 9 ,l/m b ), (7) 

with C 7 being the Wilson coefficient of the operator Q 7 . The complete set of operators 
and formulas for the Wilson coefficients can be found in Ref. [24]. a 7R , on the other hand, 
are of order l/m b [8, 9]. A p = V* s V pb and the operators Q 7 ^ are given by 

Q? R) = -^rn b sa^(l±K)bF^; 

J ( R) 

Q 7 generates left- (right-) handed photons in the decay b — > sj. The matrix element in 
(6) can be expressed in terms of the form factor Tf^ K * as 

(K*(p,r ] ) lL(R) ( q ,e)\Q^ R) \B) = 

= -^- 2 m bTl B - K \0) [^e^iiMW - (e* P )(v*q)}} 

; m b T^ K \0)S L{R) , (8) 
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where Sl,r are the helicity amplitudes corresponding to left- and right-handed photons, 
respectively, and e^(r^) is the polarisation four-vector of the photon (K*). The definition 
of Tf^ K * can be found in Ref. [25], and an updated value in Ref. [26]. In the Wolfenstein 
parametrisation of the CKM matrix, \ u ~ A 4 and is doubly Cabibbo suppressed with 
respect to A c ~ A 2 , so we drop this contribution from now on. With A M set to zero, the 
direct CP asymmetry C in (5) vanishes. 



3 Calculation of a^ R in the SM 

One contribution to A R is very well known and comes from the m s dependent part of the 
full electromagnetic dipole operator Q 7 , 

Q 7 = — [m b sa^{\ + l5 )b + m s sa^{\ - 75 )6] = Q L 7 + — Q? . (9) 
o7r m b 

Hence, a 7R is given by 

a c 7R = ^C 7 + o(^-,^-). (10) 
m b \m b m b J 

As discussed above, all contributions to A R must include a helicity flip of the s quark, 
which in the above is done by including the effects from a non-vanishing s quark mass. 



1 The a 7 ' u calculated in Ref. [20], to leading order in 1/mt, coincide with our a^ L , w±±c±ca,3 u, 7R 
zero in [20]. Our expression (6) is purely formal and does not imply that a^'^ L ^ factorise at order l/mb- 
As a matter of fact, they don't. 
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Another possibility to relieve the helicity suppression of right-handed photons is by consid- 
ering, at parton level, a three-particle final state with an additional gluon. The dominant 
contribution (with the largest Wilson coefficient) to this process comes from the operator 

Q c 2 = [^(1-75)^(1 -7 5 )&] 

and is shown in Fig. 1. As the c quark has sufficiently large virtuality in the loop (the 
photon is on-shell and the gluon nearly so), the diagram is dominated by short distances 
and can be expanded in inverse powers of m c . To do so, we follow Ref. [22] and rewrite Q c 2 

as 

Q2 = \ - 7 5 )cp7 M (l - 75)6] + 2[*7m(1 - 75) y c][crf(l - 75) y b] . (11) 

Confirming the result of Ref. [22], we find that the short-distance expansion of the diagram 
in Fig. 1 yields 

Q F = te* fM j d 4 xe l<lx T{[c(x)^c(x)}Q c 2 (0)} 

= -I^^^)^ 1 - 75)^y b} + ... (12) 

where F afS = i(q a e* /3 — q^e* ) corresponds to an outgoing photon and the dots denote terms 
of higher order in l/m c . Note that the contribution of the first term in (11) vanishes for 
an on-shell photon. The contribution of Qf to the decay amplitude is 

Aq f (B - K*>y) = -y {K^\Q F \B) , 

where 2/3 is the electric charge of the c quark and the minus sign comes from the EM 
interaction operator. At this point we would also like to make explicit our conventions for 
the strong and electromagnetic couplings. We use the covariant derivative 

D^ = d^ + ieQ fB^ - igA* y 

for a fermion with electric charge Qf. Here e = +y/Aira which is consistent with the 
sign-convention for Q 7 , Eq. (9). 2 The contribution of Q\ to a 1R is hence governed by the 
matrix element (K*-f\{DPF a ^)[s-f p {l - l 5 )gG a a ^ b]\B) , which, again following Ref. [22], 
can be parametrised as 

(K* (p, 17)7(9, e) I ( DW#) [57,(1 - 75)^y b] \B(p + q)) = 
= 2{K%p, V )\s^(l-K)9Gapb\B(p + q ))e* a qP 
= 2 {Le^e^pOq* + iL[(e*r}*)(pq) - (e*p)(r]* q )}} 

= (L + L)S L + (L-L)S R , (13) 

2 The sign of the strong coupling g differs with respect to Ref. [20], which however does not matter as 
all final expressions contain only factors g 2 . 
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where Sl,r are the photon helicity structures defined in (8). The operator Q\ thus induces 
power corrections of type (L ± L)/(mlrrib) to a c 1L and a% R , respectively. As already men- 
tioned before, these power corrections have previously been considered in Ref. [22]. Before 
we present a new calculation of L and L, let us finally give their contribution to a^ R : 

a " R = Cl ^ h ~ C2 36m c vA^(0) ' (14) 

Corrections to this expression are of order a s /mb, come with smaller (penguin) Wilson 
coefficients or are of higher order in l/mb iC . What about the convergence of the l/m c 
expansion? For the inclusive decay b — > 57 this question was studied in Ref. [27]. Higher 
terms in the short-distance expansion of (12) generate operators with higher order deriva- 
tives acting on F alS , generating powers of the photon momentum q, and on G a/ 3, generating 
new hadronic matrix elements. A complete calculation of these additional contributions to 

(14) is not possible with the presently available methods, but we can try to give an estimate. 
As found in Ref. [27], the expansion parameter of the l/m c expansion is t — (q ■ D)/(2m 2 c ) 
with D acting on the gluon field strength tensor. The hadronic matrix elements with 
additional powers of D can be estimated as 

(K*\sD n Gb\B) ~ (A QCD )" (K*\sGb\B) , (15) 

and hence t ~ (m B /2)Aq CD /(2m^) 0.2. Using (15), the l/m c series can be resummed 
and enhances the term in (12) by a factor 1.1 for t = 0.2, and 1.3 for t = 0.4. 3 Although 

(15) is only a crude estimate of the true value of these matrix elements, this result suggests 
that the l/m c expansion converges well. We also would like to mention, as noted in 
[27], that besides the derivative expansion in the gluon field there are further higher-twist 
contributions from e.g. two gluon fields. These contributions, however, are suppressed by 
additional powers of AQ CD /(m^) [27]. We shall include the effect of truncating the l/m c 
expansion by doubling the theoretical error of our final result for the CP asymmetry. 



4 Non-factorisable Soft Gluon Effects: L and L 

The following results for L and L were obtained in 1997, in Ref. [22], using the method of 
local QCD sum rules and neglecting the effects of SU (3) breaking: 

L = (0.55 ± 0.1) GeV 3 , L = (0.70 ± 0.1) GeV 3 . (16) 

Since then, a number of studies [28] have demonstrated that the appropriate method to 
calculate B decay form factors from QCD sum rules is to use QCD sum rules on the 
light-cone [29, 30]. In this paper, we cannot give any account of the method itself, but 
refer to the relevant literature. Suffice it to say that one of the main ingredients in the 
method are light-cone distribution amplitudes (DAs) of two- and three-particle Fock states 
of the final-state meson. These have been known for some time for p mesons [31, 32], but 

3 The enhancement factor is given by the function F(t,t) defined in the last reference of [27]. 
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complete expressions for K* mesons will become available only later in 2006 [33] (see also 
[34, 35]). In this paper, we include the first (preliminary) results of this ongoing study in 
our calculation of L and L. The light-cone sum rules read: 



l^B Le -ml/M* = m 4 / da 2 e~ m ^ M ^ 
m h Jo 



x 



2a\ j 



l—ao 



d(X\ 
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-,-L(2) 
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(18) 



The above expressions are accurate up to terms of order (m^./m;,) 3 , which are of higher 
twist, and 0(a s ) corrections. Here $(o) = $(0:1,0:2) are three-particle DAs of the K* 
of twist 3 or 4 (as indicated by the index). The (rather lengthy) definition of these DAs 
is given in Ref. [33]. The variable a± can be interpreted as the longitudinal momentum 
fraction carried by the quark in the meson, whereas a 2 is the momentum fraction carried 
by the antiquark. I [$] is a functional acting on the DA $(o), which is defined as 



a 2 



1-x 



dx / da\^{a\, x). 
'0 Jo 

The DAs can be described in a systematic way using conformal expansion [31]; here, 
we restrict ourselves to the leading terms in that expansion and use the expressions [33] 
(o 3 = 1 — a 1 — a 2 ) 

$!U*(«) = 360oio 2 a^ + ul K (an - a 2 ) + (7a 3 - 3)| , 



%,K*(®) = 360oia 2 «3 |Csk + KA a i - "2) + £3^2 ( 7as ~ 3 ) J ' 
$i;K*(QL) = 360a 1 a 2 o 3 J j/c^- + ^(«i - a 2 ) + \j K ^ (7o 3 - 3)| , 

(a) = 120OiO 2 «3 Q «3A- + ^ «4k) , 
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H = 1 GeV 


A 

^3K 


0.033 ±0.007 


^3K 


0.06 ±0.03 


~ll 
UJ 3K 


-0.06 ±0.02 


K 3K 


0.001 ±0.001 


U 3K 


0.14 ±0.03 


^3K 


-0.02 ±0.01 


K 3K 


0.006 ±0.003 


U 3K 


0.4 ±0.1 


^3K 


-0.05 ±0.02 


Cik 


0.10 ±0.05 




= — Cik 


K 4K 


0.012 ±0.004 



Remarks 



new; Q p determined in [37] 
G-odd, new 

determined in [37] 



new; to 



3p 

G-odd; previously determined in [38] 
new; u\ p determined in [37] 
G-odd, new 
G-odd, new 

new; determined in [31] 
G-odd, new 

quoted from [29]; no SU(3) breaking; to be updated in [33] 
quoted from [29]; no SU(3) breaking; to be updated in [33] 
G-odd; quoted from [35] 

Table 1: Three-particle twist-3 and 4 hadronic parameters. All results labelled "new" are pre- 
liminary and will be finalised in Ref. [33]. Note that the absolute sign of all these parameters 
depends on the sign convention chosen for the strong coupling g. The above results correspond 
to the choice D M = <9 M — igA p (\ a /2) of the covariant derivative. 
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1 a 1 



(ai - a 2 )(i K > , 



(19) 



Preliminary numerical results for the various hadronic parameters (, k, uj and A are col- 
lected in Tab. 1; they will be discussed in more detail in Ref. [33]. The DAs defined above 
are related to those introduced in Ref. [31, 32] as 



(20) 
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Although the introduction of new notations may, at first, look unmotivated, it actually 
extends the labelling scheme introduced, in Ref. [36], for pseudoscalar mesons, to vec- 
tor mesons and aims to provide a systematic way to label the multitude of two- and 
three-particle pseudoscalar and vector meson DAs, replacing the slightly ad-hoc notations 
introduced in our previous papers on the subject [31, 32]. As for the other hadronic pa- 
rameters entering (17) and (18), we use m b = (4.7 ± 0.1) GeV, f B = (200 ± 30) MeV, 
fl, = (217 ±5) MeV [39] and /£.(lGeV) = (185 ± 10) MeV [34]. All scale-dependent 
parameters are evaluated at the scale /i 2 = m% — m\ ± 1 GeV 2 , see Ref. [25]. s and M 2 
are sum rule specific parameters which do not acquire sharp values, but have to be varied 
in a certain range. Based on our experience with B decay form factors [25] we choose 
s = (35 ± 2) GeV 2 and M 2 = (10 ± 3) GeV 2 . We then obtain 

L = (0.2 ± 0.1) GeV 3 , L = (0.3 ± 0.2) GeV 3 , L — L — -(0.1 ± 0.1) GeV 3 . (21) 

It turns out that the contribution of the (m K */m b ) 2 terms to the sum rules is tiny, so that 
the result and its uncertainty is entirely dominated by m b , fs and the twist-3 DAs §l. K * 
and We repeat that the parameters describing these DAs, collected in Tab. 1, are 

preliminary. 

The results in (21) refer to the renormalisation scale fi 2 = m 2 B — m\ m (2.2 GeV) 2 . 
Unfortunately, the dependence of L and L on /j, is unknown. We can, however, estimate 
the potential impact of a change of scale by evaluating the light-cone sum rules at the higher 
scale ji = m b , although this is, strictly speaking, incorrect in that framework. Nonetheless, 
L and L itself decrease by about 20% by this procedure, whereas L — L decreases by 10%, 
which is well within the quoted errors. 

Comparing with the results obtained in Ref. [22], Eq. (16), we find that our central 
values are considerably smaller. As mentioned before, the authors of [22] used local sum 
rules, which are of only limited value for determining B decay form factors at maximum 
recoil, i.e. for maximum energy of the final state meson, see the discussion in the first two 
references in [28]. On the other hand, in 1997 not much was known about three-particle 
twist-3 DAs of vector mesons, so local QCD sum rules were the best tool at hand at the 
time. We also find that our errors are larger than those in (16), which is due to the fact 
that the uncertainties quoted in (16) are obtained by varying only the sum rule parameters 
s and M 2 , but not the hadronic input parameters. 

5 Results and Conclusions 

We are now finally ready to present results for the CP asymmetry S in (1). The m s - 
dependent terms in (3) yield 

S SM ' SR = -0.027 ± 0.006(m Sj6 ) ± 0.001(sin(2/3)) , (22) 

where we use m s (2GeV) = (100 ± 20) MeV [40], m b (m b ) = (4.20 ± 0.04) GeV [41] and 
sin(2/3) = 0.685 ± 0.032 [16]. One can estimate the impact of radiative corrections on that 
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result by comparing it with the perturbative QCD calculation of Ref. [21]. The authors of 
[21] obtain the same central value for S SM ' SR and also quote, very helpfully, results obtained 
for neglecting various sources of corrections, in particular — 0.034±0.013 if all long-distance 
contributions are neglected. From this we conclude that the impact of radiative corrections 
on (22) is likely to slightly increase the asymmetry, but not by more than 0.01. As for the 
contribution of L — L, it is given, to leading order in a s , by 

s sm,„„ glm , = _ 2sin(2/J) J_| = . 005 ± . M . (23) 

Here we use C2(mb) = 1.02, CV(m{,) = —0.31, which are the leading-order values, m c = 
1.3 GeV and Tf^ K *(0) = 0.31 ± 0.04 [26], and have doubled the error to account for 
neglected higher-order terms in the l/m c expansion. That is: the contribution of soft 
gluons to S is much smaller numerically than that in m s /mb, Eq. (22). This result has 
to be compared with the dimensional estimate presented in Ref. [9], from a SCET-based 
analysis, 

|£SM,soft g iuons| = 2sin(2/3) ( 



3CV 



A QCD _ Q6 _ (24) 



Our result (23) suggests that the true value of the soft gluon contributions is much smaller. 
Comparing (23) and (24), it becomes obvious that this is mainly due to the factor 1/36 in 
(23) resulting from the short- distance expansion of the charm loop in Fig. 1. 

While the aim of our letter was to calculate the soft gluon contributions to S SM induced 
by the operator Q 2 and to check the estimate of Ref. [9] that it could induce a 10% effect, 
our result for S SM ' soft s luons has now become, due to the suppression factor 1/36, that 
small that one may start to wonder about the size of other corrections. One source of 
such corrections, and actually probably the dominant one, are radiative corrections to the 
term in m s /mb which we estimate using the results of Ref. [21] in the way discussed above. 
Another class of (soft gluon) corrections are diagrams with the same topology as Fig. 1, but 
a different operator. As long as there is a charm quark in the loop, these contributions are 
controlled by the matrix element L—L, but suppressed by small penguin Wilson-coefficients 
Cpeng < 0.1 and hence can be neglected. For light quarks in the loop, one cannot apply 
the short- distance expansion as done in this letter, but has to follow a different approach. 
We will discuss this approach in a forthcoming paper on power-corrections to B — > fry 
[42]; the result is that the contribution of light quark loops is of approximately the same 
size as that of charm loops, so again these contributions are suppressed by small Wilson 
coefficients. A second, different topology is given by annihilation diagrams induced by the 
penguin operators (sb)v-A{dd)v±A- This contribution is enhanced by the fact that it is 
a tree diagram; it can be calculated using the results obtained in Ref. [43] for B — > •yev 
transitions. For the contribution with the largest Wilson-coefficient from the penguin 
operator Q^, one has 

a m - a 7R + C 4 - m s T B^. (0) (*V(0) " F A (0)) , 
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where Q u ^ are the electric charges of the corresponding quarks and Fy,A are the form factors 
determining the B -> 7 transition. Using F v (0) - F A (0) « 0.016 [43] and C 4 (m b ) = 0.08, 
the shift of a^ R turns out to be ~ 0.3 • 10~ 3 which is to be compared with the (dominant) 
m s /mb term m 6 • 10~ 3 and the term in C2: ~ 1 • 10~ 3 . Let us note in passing that 
Fy{0) — Fa{0) is induced by long-distance photon emission and given in terms of three- 
particle Fock states of the photon [44], so also for this contribution the necessary spin-flip 
in the parton-level process b — > S7 is induced by a higher Fock state, this time of the 
photon. One more possible topology are hard-spectator scattering diagrams involving the 
chromomagnetic dipole operator Q$. Although we cannot give a firm estimate of this 
contribution to a^ R , we expect it to come mainly from long-distance photon emission 
governed by the same three-particle Fock state of the photon mentioned before and to 
contribute at the same level as the other terms discussed above. Although we were not 
able to identify a further sizable contribution we will add an uncertainty of 0.01 to the 
asymmetry. 

Our final result for the CP asymmetry is the sum of (22) and (23) with uncertainties 
added in quadrature: 

gSM = S SM,s R + ^SM.soft gluons = _ . 22 ± 0.015+° 01 , (25) 

where the second uncertainty accounts for neglected contributions induced by penguin 
operators and the chromomagnetic dipole operator, and the third, asymmetric uncertainty 
is to account for neglected 0(a s ) corrections to (22), which, based on the results of Ref. [21], 
we estimate to be negative and not to exceed 0.01. In principle these corrections can also 
be calculated in QCD factorisation, but this goes beyond the scope of this letter. 

To summarize, we have calculated the dominant contributions to the SM prediction 
for the time- dependent CP asymmetry S in B° — > K* ^. These come, on the one hand, 
from terms in m s /m b , and on the other hand from short- distance processes involving an 
additional gluon, see Fig. 1. We find that in contrast to recent suggestions that the latter 
be large, they are actually substantially smaller than the former. Additional hadronic 
corrections to our result are expected to be even smaller and due to radiative corrections, 
small Wilson coefficients and higher order terms in the heavy quark expansion. The most 
dominant correction is likely to be radiative corrections to (22), which have already been 
calculated in perturbative QCD and found to be —0.01. A confirmation of this result in 
QCD factorisation, if possible, would be welcome. Our result (25) confirms that the CP 
asymmetry is an excellent quasi null test of the SM in the sense of Ref. [11] and that any 
significant deviation of the experimental result from zero will provide a clean signal for 
new physics. 
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